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lO . Abstract 

o 

' Starting from the QCD Hamiltonian in near-light cone coordinates, we study the 

■ dynamics of the gluonic zero modes. EucUdean 2+1 dimensional lattice simulations 

, show that the gap at strong coupling vanishes at intermediate coupling. This result 

I opens the possibility to synchronize the continuum limit with the approach to the light 



cone. 



5^ ■ 1 Introduction 



The solution of Quantum Chromodynamics on the light cone is still an unsolved theoretical 
^ ■ task the present status of which is reviewed in ref. In a recent paper [Q] a formulation 
■ of QCD in coordinates near the light cone has been proposed which has the advantage 



of keeping a direct link to equal time theories. The Cauchy problem is well defined 
in near light cone coordinates, since the initial data are given on a space like surface. This 
formulation avoids the solution of constraint equations which, on the quantum field theoretic 
level, may be very complicated. The problem of a nontrivial vacuum appears in a solvable 
form related to the transverse dynamics. This is physically very appealing, since in high 
energy reactions the incoming particles propagate near the light cone and interact mainly 
exchanging particles with transverse momenta. We would like to connect successful models 
for the soft transverse nonperturbative dynamics of high energy reactions (c/. refs. [§,0]) to 
the underlying QCD Hamiltonian. In diffractive reactions, fast Lorentz contracted hadrons 
experience the confining forces of QCD in their transverse extensions and interact with 
other fast moving hadrons via soft interactions. In a theory of total cross sections the 
nonperturbative infrared dynamics in the transverse plane is essential. The objective of this 
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paper is to investigate the effective transverse Hamiltonian in QCD near the hght cone for 
small light cone momenta. 

Near light cone QCD has a nontrivial vacuum which we claim cannot be neglected even 
in the light cone limit. We will demonstrate the existence of massless excitations in the zero 
mode theory. These excitations do not decouple in the light cone limit. Genuine nonper- 
turbative techniques must be used to investigate the behavior of this limit. In principle the 
additional parameter which labels the coordinate system can be chosen arbitrarily. We will 
show that the zero mode Hamiltonian depends on an effective coupling constant containing 
this parameter and evolves towards an infrared fixed point. Therefore, we propose a sophis- 
ticated choice of the frame dependence which facilitates calculations in the near light cone 
frame. We use the infrared fixed point of the zero mode Hamiltonian to follow a trajectory 
in the space of couplings, where high resolution is synchronized with the light cone limit. 
Note, we consider the parameter associated with the frame dependence as yet another cou- 
pling constant. As shown in Ref. the zero mode sector is also relevant to the definition 
of M-theory in light cone coordinates. 

We choose the following near light cone coordinates which smoothly interpolate between 
the Lorentz and light front coordinates: 




x" = ^ X -x^ . 1 



The transverse coordinates x^, x^ are unchanged; x* = x^ is the new time coordinate, x~ is 
a spatial coordinate. As finite quantization volume we will take a torus and its extension in 
as well as in "1, 2" direction is L. The scalar product of two 4- vectors x and y is given 
with x_Ly_L = x^y^ + x^y'^ as 

^iiV^ = x^y^ + x^y^ — rfx^y^ — x±y± 

= x-y+ + x+y- + rfx+y+ - x±y± . (2) 

Obviously, the light-cone is approached as the parameter t] goes to zero. For non-zero rj, the 
transition to the coordinates introduced above can be formally identified as a Lorentz-boost 
combined with a linear transformation, which avoids time dependent boundary conditions [Q. 
The boost parameter /5 = f 3 is given by 

1 -|- ?7^/2 

indicating that for 77^ — »• the relative velocity V3 —>■ c{= 1). Of course, this is connected 
to the well-known interpretation of the 'tilted' light-cone frame in terms of the infinite 
momentum frame. 

The choice of near light cone coordinates allows to quantize the theory on a space-like 
finite interval of length L in x~ at equal times, i.e. Ax'^ = 0. The invariant length squared 
of this interval for Ax^ = is related to the length of the compact x~ dimension. 

As^ = A,T"Ax+ + A,T+A,T" - r7^(Ax")^ - Ax^ 

= -V'L\ (4) 
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For simplicity we consider also the transverse dimensions periodic in L. Previous work of 
the St. Petersburg and Erlangen groups assumed a fixed tilted coordinate system with 
fixed transverse ultraviolet cut off. Our main purpose in this paper is to consider the zero 
mode fields on a transverse lattice with varying transverse lattice spacing a. We propose 
to approach light-cone dynamics by synchronizing the continuum limit A = vr/a — > oo with 
the light-cone limit ^(^° ^^)' "^^^ vacuum fiuctuations in the transverse directions 

induce a second order phase transition which allows to have a continuum limit of the lattice 
theory. Thereby we can eliminate the cut off in a controlled way, preserving the nontrivial 
vacuum structure. 

The gauge fixing procedure in the modified light-cone gauge d-A_ = involves zero 
modes dependent on the transverse coordinates. These zero mode fields carry zero linear 
momentum p_ in near light cone coordinates, but finite amount of po + Ps- They correspond 
to "wee" partons in the language of the original parton model of Feynman. In SU{2) the zero 
mode fields a_{x±) are proportional to r^, i.e. they can be chosen color diagonal. The use of 
an axial gauge is very natural for the light-cone Hamiltonian even more so than in the equal- 
time Hamiltonian. The asymmetry of the background zero mode naturally coincides with the 
asymmetry of the space coordinates on the light cone. The zero mode fields describe disorder 
fields. Depending on the effective coupling the zero mode transverse system will be in the 
massive or massless phase. Its second order phase transition allows to perform the continuum 
limit in the zero mode Hamiltonian. Our main conjecture is that the continuum limit and 
light cone limit can be realized simultaneously at this critical point. The evolution of the 
coupling determines the approach to the light cone with transverse resolution approaching 
zero. The resulting relation is reminiscent of the simple behavior one gets from considering 
naive scaling relations and dynamics in the infinite momentum frame. We will also calculate 
the contribution of the two dimensional zero modes to the ground state energy and show 
that it scales with the three dimensional volume. This estimate demonstrates very simply 
the relevance of modes with lower dimensionality to the full problem in the case of critical 
second order behavior. 

The naive considerations for a simultaneous light cone and continuum limit go as follows. 
We characterize the 'infinite' momentum frame by giving the momenta of the proton and 
photon in usual coordinates. The fast moving proton carries {P,0±,P) and the photon 
q = VQ^, —^), where u is the energy transfer in the laboratory. The dominant 

contribution comes from energy conserving transitions, where the energy of the quark in the 
final state equals the sum of initial quark and photon energies. For collinear quarks with 
momentum p = {xbP,0±,xbP), where xb denotes the momentum fraction, one explicitly 
finds 

\IQ'+[^bP-^) =^bP + ^, (5) 

which yields the "scaling variable" xb = as long as P > Q/xb, i-e. the "infinite" 
momentum has to be large enough: 

-fM = P >Q = 1/a. (6) 

At the same time, the photon has a wavelength (x 1/Q which resolves transverse details of 
size a in the hadronic wavefunction. By eq. (^, the 7 factor of the "infinite" momentum is 
related to rj. Therefore we expect that the parameter 77 approaches zero with the transverse 
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cutoff A = vr/a — s> oo 

One of the objectives of this paper is to derive the precise relation between 77 and a. We will 
simulate the zero mode dynamics on a (2 + 1) dimensional lattice to find the fixed point of 
the effective coupling. The light-like limit is governed by an infrared fixed point (c/. ref. §). 
Most previous discussions of the zero mode problem have been on the classical equations of 
motion level, here we will work on the fully nonperturbative quantum level. Our work is not 
directly related to M-theory, but some of the consequences may be relevant for other studies 
with light like compactification. 



2 Near Light Cone QCD Hamiltonian 

In Ref. the near light cone Hamiltonian has been derived. A similar Hamiltonian has been 
obtained in equal time coordinates [^ll^]- Here we will sketch the derivation. We restrict 
ourselves to the color gauge group SU (2) and dynamical gluons; only an external (fermionic) 
charge density pm is considered here. 

Since the coordinates have no momenta conjugate to them, the Weyl gauge A*]^ = 
is the starting point for a canonical formulation. The canonical momenta of the dynamical 
fields A"i , A^ are given by 

djC 

n- = -Q^ = F%^^^Fl, (8) 
^From this, we get the Weyl gauge Hamiltonian density 

Uw = ^n«_n"_ + + E (n- - f^,)' . (9) 

The Hamiltonian has to be supplemented by the original Euler-Lagrange equation for A^ 
as constraints on the physical states (Gauss' Law constraints) 

G\x^,x-m = {Dl'ut+Dfu'^ + gpl)\<l>) 

' D":^Ut + Gl) \<^) = . (10) 



In order to obtain a Hamiltonian formulated in terms of unconstrained variables, one 
has to resolve the Gauss' Law constraint. Via unitary gauge fixing transformations [9,10] a 
solution of Gauss' Law with respect to components of the chromo-electric field H_ can be 
accomplished. This gives a Hamiltonian independent of the conjugate gauge fields A^, i.e. 
the latter become cyclic variables. Classically this would correspond to the light front gauge 
y4_ = 0. However, this choice is not legitimate if we want to consider the theory in a finite 
box. Instead, the (classical) Coulomb light front gauge d-A^ = is compatible with gauge 
invariance and periodic boundary conditions. The reason is that A_ carries information on 
the (gauge invariant) eigenvalues of the spatial Polyakov line matrix 



V{x±) = Pexp 



ig J dx A_{x±,x 



(11) 
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which can be written in terms of a diagonal matrix with the zero mode field a^{x±)Y = 

V[xa_) = V exp [igLa. {x±)] . (12) 

Obviously we have to keep these 'zero modes' a^{x±) as dynamical variables, while the other 
components of A_ are eliminated. The zero mode degrees of freedom are independent of x~ 
and, therefore, correspond to quantities with zero longitudinal momentum p_ . 

In order to eliminate the momentum n_, conjugate to A_, by means of Gauss' Law, one 
needs to 'invert' the covariant derivative D_ which simplifies to d- = d- — ig[a-, ..]. On the 
space of physical states one can simply make the replacement]^ 

n_(x_L,x_) p-ix±) - (dZ^) G_i_{x_i_,y~) . (13) 

The operator p-{x±) is also diagonal and p^{x±) is the momentum conjugate to the zero 
mode a^{x±). It has eigenvalue zero with respect to d_, i.e. d_p_ = 0, and is therefore not 
constrained. 

The appearance of the zero modes implies a residual Gauss' Law, which arises from the 
x_ integration over eq. (p!OD for the a = 3 component using periodic boundary conditions in 
the X- direction. This constraint on two dimensional fields can be handled in full analogy 
to QED, since it concerns the diagonal part of color space. A further Coulomb gauge fixing 
in the SU{2) 3-direction eliminates the color neutral, ^"-independent, two-dimensional 
longitudinal gauge field in favor of a neutral chromo-electric field 

3 

e±ix^) = g\/± I dy~dy^d{x^-y^) [f^' Al{y^,y~)Ii\{y^,y-) + pl{y^,y-)] y . (14) 
Here we use the periodic Greens function of the two dimensional Laplace operator 

d{z^) = E ie-"- , p. = |n , (15) 

where n = (ni,n2) and ni,n2 are integers. 

As a remnant of the local Gauss' Law constraints, a global condition 

Q'\^') = J dy-dy^ {f^''Al{y^,y-)Ui{y^,y-)+pl{y^,y-)} |$') = (16) 

emerges. The physical meaning of this equation is that the neutral component of the total 
color charge, including external matter as well as gluonic contributions, must vanish in the 
sector of physical states. It can be shown (see 10) that there must be Qi2{x±) = everywhere 
in transverse space in order to avoid an infinite Coulomb energy. The two conditions together 
suggest that physical states have to be color singlets. 

The final Hamiltonian density in the physical sector explicitly reads (in terms of uncon- 
strained A± and n_L obtained after a shift is made by subtracting averages) 

n = tT[diA2-d2Ai-tg[AuA2f + \tT[U^-{d-A^-ig[a^,A^])f 



1 ^2 



+ ^tr 



+ ^P^-\x±)pi{x±) 



+ ^ f\,- f\y- ^ ' G^,p{x^,Z-)G^,,ix^,y-) ^ ^,2.n(.-^y-)/L ^ ^^j^ 



^° [^ + g{a^g{x^)-a^,{x^)) 



^The inversion of d- can be explicitly performed in terms of its eigenfunctions, cf. (pO|). 
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where p and q are matrix labels for rows and columns, a_g = {cL-)qq and the prime indicates 
that the summation is restricted to n 7^ if p = g. The operator G± {x±,x~) is defined as 

= V^n^ + gr'^Y"^", (ni - iel) + gpm • (18) 

The last two terms of the Weyl gauge Hamiltonian come from the original term 11^ with 
squared electric field strengths in x~ direction. After elimination of n_ the zero mode part 
and the light cone Coulomb energy in the axial gauge remain. In the Coulomb term one sees 
the role of the zero mode fields as infrared regulators of the spatial momenta p_ = 27m/ L 
which are quantized due to the compact interval L. Since the two dimensional theory has 
been Coulomb gauge fixed, the electric field e± replaces the canonical momentum of the 
longitudinal, neutral gauge field which has been eliminated. 

We note that the terms containing and the momentum n^, have the pre-factor l/rj'^. 
Physically this pre-factor signals the increase of transverse electric energies with the boost 
factor 7 = (\^r])~^. The boost also couples transverse electric fields with transverse magnetic 
fields. In the light-cone limit the pre-factor diverges and the adjacent brackets become 
constraint equations. This reflects the corresponding reduction of the number of degrees of 
freedom if one goes exactly on the light cone. We do not follow this procedure, but keep a 
finite 77 as a kind of Lagrange parameter. 

A characteristic feature of an exact light-cone formulation is the triviality of the ground 
state. This may simplify explicit calculations, e.g. of the hadron spectrum. However, the 
light cone vacuum is definitely not trivial in the zero mode sector. In fact, already in ref 
we have shown strong and weak coupling solutions of the zero mode Hamiltonian. Massless 
modes influence the dynamics in the light cone limit, whereas massive modes decouple when 
rj ^ 0. In contrast to earlier work, we solve the zero mode Hamiltonian in this paper 
numerically showing the transition from the massive phase to the massless phase. 

The zero mode degrees of freedom a_ couple to the transverse three dimensional gluon 
fields Ai via the second magnetic term in Ti, the Coulomb term and directly via the electric 
field e±. We remark that the transverse electric fields H_|_ and e± are dual to the magnetic 
fields d_A± and V_La_. This duality is typical for the light cone and is absent in the equal 
time case. Since duality plays an important role in supersymmetric QCD its role in light 
cone theories should deserves to be investigated in greater detail. 

In the following we neglect the couplings between the three and two dimensional fields 
and consider the pure zero mode Hamiltonian 



h= d^x 



(19) 



The global constraint, eq. ([T^), does not contain a_ and p_ and, consequently, is irrelevant 
for the time being. Even at this level of severe approximations the zero mode Hamiltonian 
differs from the corresponding one in QED. The reason is the hermiticity defect of the 
canonical momentum p~ . In the Schrodinger representation eq. (|19D 

J ~2LJ(aUx^ )) daUx^ ) T-^^^^^J dalix^ ) ^ 2^^^ ' 



h 



(20) 



2LJ{a^_{x^))daUx^) ^ ^'^ dal{x^) 2r]^ 
contains the Jacobian J{a_) which equals the Haar measure of SU(2) 

J (al{x^)) = sin^ ( ^alix^)) . (21) 
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It stems from the gauge fixing procedure, effectively introducing curvilinear coordinates. It 
also appears in the functional integration volume element for calculating matrix elements. 
It is convenient to introduce dimensionless variables 

gLal{x±) 

(p{x±) = , (22) 

which vary in a compact domain < (/? < tt. We regularize the above Hamiltonian h by 
introducing a lattice spacing a between transversal lattice points b. Next we appeal to the 
physics of the infinite momentum frame and factorize the reduced true energy from the 
Lorentz boost factor 7 = \/2/ rj and the cut off by defining h^ed 

h = -^/ired (23) 

2r]a 

In the continuum limit of the transverse lattice theory we let a go to zero. For small lattice 
spacing we obtain the reduced Hamiltonian 

-^j:{^ib)-fib + e.)f] . (24) 

9eS 1=1,2 ) 

.i^^ (25) 

The first part of the Hamiltonian contains the kinetic (electric) energy of the SU{2) rotators 
on a half circle at each lattice point and the second part gives the potential (magnetic) 
energy of these rotators due to the differences of angles at nearest neighbor sites. For further 
discussion we define these electric and magnetic parts as: 

/ired = K + hra (26) 
he = E^l 



^red ~ 

with the effective coupling constant 




i=l,2 



byb+Si 



In the strong coupling domain where we have analytical solutions of the zero mode Hamil- 
tonian, the numerical Hamiltonian lattice theory agrees with the analytical solutions for the 
mass gap. The real question concerns the continuum limit of the zero mode system which 
occurs outside of the strong coupling region. For this purpose we have to find the region of 
vanishing mass gap for the lattice Hamiltonian. 



3 Lattice Calculation of the Zero Mode Hamiltonian 

We solve the equivalent lattice theory in a Euchdean formulation. The zero mode Hamilto- 
nian represents a 2 -|- 1 dimensional theory in two spatial and one time direction. The lattice 

has N^a = Nya extensions in transverse space and Nt^t = T extension in near light cone 
time. To set up the density matrix one has to write down the Trotter formula for the given 
Hamiltonian. Using h,:cd = he + h^n wc have 

exp(-T/ired) = ,lim [exp {-AThj2) exp {-Arh,) exp {-AThj2)f^ , (27) 
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where each time evolution step Ar can be separately done for the electric and magnetic 
part of the Hamiltonian. For definiteness, we choose Ar = a/2 in all the following. In 
the Appendix we show that this choice of Ar optimizes the updating procedure since it 
generates approximately equal widths for the weight functions resulting from the kinetic 
and potential energies. The different time slices will labeled with the index I. The electric 
Hamiltonian can be evaluated by inserting products of complete set single site eigenfunctions 
Cm {cf. 0). Practically a maximal number iVmax = 100 of eigenfunctions is fully sufficient 
to reach convergence in the interval of couplings we need. 

ax 

ni=0 

with the (single site) eigenfunctions and eigenvalues given as: 




CnM) = g'eS niim + 2) CnM) ■ (30) 



The eigenfunctions form an orthonormal set with respect to a scalar product which contains 
the Jacobian in the measure. The Jacobian J (^ip(b)^ has been defined above (pTj). The 

magnetic part is diagonal in {(f{b)} i.e. local in time: 

iWimhr^iWim = huwim = E E ^r^+e. H^) - ^.{Ue, )) . (si) 

The full partition function is given by an integral over all time slices 

Z = trexp[-T/ired] =/n(^(<^W)M^)) ({v^W}|exp[-TM|{v^(&)}) (32) 

b 

(A^max 
J(^(6))rf^(6) Cni^i+i{b))exp[-gl^ni{ni + 2)AT]CnMb)) 
ni=0 

X l[exp[-h^{{ipi{b)})AT] . 
I 

Because of the Jacobian the dominant contributions to the partition function come from 
values around |. The Hamiltonian is invariant under refiections yj; — | ^ ^ — In the 
strong coupling limit, where (7^^ is large, the half rotators act almost independently on each 
lattice site producing a large mass gap. The system is disordered. For decreasing coupling 
constant g'^g the movement of the individual rotators becomes locked from one site to the 
next. Long range correlations develop. The order parameter (or "magnetization") of the 
system is the expectation value of the trace of the Polyakov line (|TT|), which on the lattice 
has the form 

P=ltrV = ^Y^os^{b,T) . (33) 

b,T 

The order parameter is odd under the above symmetry operation. Operator expectation 
values are evaluated with the density matrix defined by h^ed and T: 

(0) = ltr(Oexp[-TM) • (34) 
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At each particular (3g 

^ = l/^7efr (35) 

we simulate lattice sizes with equal number of sites in all directions = = Ny = Nt 
and = 4, 6, 8, 12, 16. A Metropolis algorithm is used for updating with the steps size and 
the number of hits adapted to (3g. In order to tabulate the Boltzmann weights related to 
the timelike links we also discretized the continuous angle variables to a system of A^^ = 60 
orientations: Lpj = jn/N^, j = [0, A";^]. 

In this explorative investigation we generated between 5000 and 50000 uncorrected con- 
figurations depending on the Pg and lattice size. All our calculations were done on a cluster 
of AlphaStations with the single site Metropolis updating algorithm sketched above, so the 
accuracy can be improved using a more powerful algorithms and/or more computing re- 
sources. 

We calculate the Pg dependence of the following quantities: average electric energy Se 
and magnetic energy per site, the average of the absolute value of the order parameter 
(|P|), the susceptibility x ^i-nd the normalized fourth cumulant Qr (which gives the deviation 
of the moments of the Polyakov expectation value from a pure Gaussian behavior): 

^e,m = -^{he,m) (36) 

X = N'{{P')-{Pr) (37) 
9r = ^-3. (38) 

Firstly we calculated the ground state energy for strong couplings g'^g and compared 
with the exact calculation 0. The obtained results from the lattice agree with the analytical 
result. For higher values of jSg the lattice calculation agrees within 10—20 % with the previous 
effective double site calculation for the energy per site 0. The ultraviolet regularization 
with At influences the final lattice result. The energies are measured with high accuracy 
and practically do not depend on the infrared cutoff, i.e. on the total lattice size for A^ > 8. 
The situation for other variables is different. 

In order to see the emergence of a massless phase, we investigated the (connected) time 
correlation functions of the (trace of the) Polyakov line operators: 

mr) = ^ (Ecos(^(^0))cos(^(&',r))) - {\P\)' . (39) 

\ b,b' ' 

The correlation masses are obtained from a fit of the time correlation functions ^'(O, 'r)|^=„Ar 
to a parameterization taking the periodicity in time into account 



K{0,nAT) = ccosh 



mAr n 

2 



(40) 



The obtained masses for different lattice sizes are shown in Fig.|l]a. For each fixed jSg we 
extrapolated the mass to the limit N —>■ oo using a linear 1/A^ parameterization 

m{(3g,N)=mo{f3g) + ^. (41) 
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Figure 1: Correlation masses m for different effective couplings f3g = l/f^efr- (^) ^ 
different lattice sizes N from a fit of the measured correlation functions K{0, r) with the 
parameterization (|40D; (b) - in the limit N ^ oo (see text). The solid line on (b) is the fit 
curve given by expression (^21) . 



Our data for ■m{(3g,N) can be fitted with this linear form ( p] ) for f3g^3. The resulting 
dependence of mo on f3g is presented in Fig.|pD. The mass vanishes in the infinite volume 
limit near /3 ~ 5, which is a first guess for the critical coupling P*. 
In order to extract more exact information on the critical behavior 



3* 



(42) 



from our data on small systems we have done a Finite Size Scaling (FSS) ||T2,13| analysis of 
the variables (|-P|), x and gr, searching for the critical coupling /5* in the interval 3 < Pg < 9. 
We could determine the critical indices /3, 7 and z/ using an approach employed for SU{2) 
gauge theory at the finite temperature transition by Engels et al. |T^. The general form of 
the scaling relations for a variable V {V = {\P\), Xi 9r) is 



(43) 



where p is the corresponding critical index (/?, 7, 0) for the respective quantities and t is the 
reduced inverse coupling [3g: 

^ (44) 



t 



In practice eq. (^31) is computed near t = 0. Expanding up to first order in t and taking into 
account only the largest irrelevant exponent yi = —uo one obtains 



y (t, N) = N^/" [co + (ci + C2N-^) tN^/" + CgiV^ 



(45) 
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Figure 2: Comparison of the theoretical parameterizations given in eqs. (44-45) for (|P|) 
and X with the lattice data results as a function of (5g for different lattice sizes A^. 



With this parameterization we analyze our data for lattice sizes = 4, 6, 8, 12, 16. The 
/3g dependence in eq. ( ^5l ) is valid near the critical point. The quality of the fit to the data 



in different intervals on /3g can serve as a guide to localize the critical point /5* . We find 
that in our case this parameterization can be applied for and x only for > 5: for 
smaller f3g the x^/D.F. becomes too large. Finally we use the parameterization of eq. ( ^Sf ) 
in the interval 5 < /?g < 9 (see Fig.^j) where |t| ^ 0.3 and an expansion of first order in t 
still is applicable. In this [5g interval the fits for and x yield a x^/D.F ~ 2. To have 
a reliable error estimate for the extracted parameters we use the jackknife method. The 
results for /?* and critical indices /?, 7 and z/ are presented in Table 1. We use a; = 1, but in 
fact the results do not depend on uo in the interval 0.8 <uj< 1.2. 

It should be noticed that the system frequently flips between the two ordered states on 
a finite lattice. Therefore the expectation value (P) in equation (|37|) vanishes with good 
accuracy. Analogously to the treatment of the magnetization in the 3-dimensional Ising 
model, we should - instead of the expression for x eq. ( p7| ) - define the susceptibility in 
the broken phase as Xbrokcn = {{P^) — (|P|)^). For the Ising case this susceptibility is 
supposed to converge towards the correct infinite volume limit. Due to lack of statistics, 
however, we did not use separate expressions for x in the two phases. Our data for gr are 
even much less accurate which excludes the possibility to use them in the FSS analysis. 

The values found for the critical indices /3, 7 and u are not far from those of the Ising 
model {D = 3) (see Table 1). They are much less accurate but also in agreement with a 
high statistics analysis ||14[ of the SU{2) deconfinement transition at finite temperature. We 
interprete the behavior of the light cone zero mode theory as a consequence of the underlying 
Z{2) symmetry of the Hamiltonian. Z{2) transformations correspond to reflections of the 
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V 




Transverse QCD (2 + 1) 


Ising {D = 3) 


^ 1 Pi \ 


o 
P 


n 01 -1- n r\R 
U.zi ± U.Uo 


U.ozo ( (^iUj 




V 


0.57 ± 0.06 


0.6289(8) 






6.3 ± 1.1 




X 


7 


i.iO IC u.oo 








0.52 ± 0.10 


0.6289(8) 




PI 


5.7 ± 1.2 





Table 1: Critical parameters for the transverse QCD near the light cone. The Ising indices 
are from ref. [03 . 



angle variables cp around ^. The Hamiltonian and the measure of eq. (17) are invariant 
under these transformations. The resulting critical behavior is common to Ising like models. 
The critical coupling itself is a nonuniversal quantity and we have found a rough value. 
More extended work with the lattice Hamiltonian near the light cone is needed to clarify the 
continuous transition further. 



4 Conclusions 

The scaling analysis gives indications that there is a second order transition between a phase 
with massive excitations at strong coupling and a phase with massless excitations in weak 
coupling. To reach a higher accuracy large scale simulations of the Hamiltonian zero mode 
system are needed. In this context it is advisable to treat the coupled system of three and two 
dimensional modes together. A calculation in the epsilon expansion |]ll| gives the zero of the 



/3-function as an infrared stable fixed point. This shows that the limit of large longitudinal 
dimensions L is well defined. Using the running coupling constant g'^g of the zero mode 
system we have (c/. (^2])) 

= (als - 9li) (46) 

Cofl'efr 

with gll = 0.17 ± 0.03 and z/ = 0.56 ± 0.05 from our lattice calculations (c/. Table The 
mass can be interpreted as an inverse correlation length ma = a/^. When the correlation 
length ^ approaches the temporal size of the lattice ^ ~ Nxa/'i in the infrared limit we can 
identify ma with a/L, where L is our infrared length scale. Therefore the effective coupling 
9eS = 5'^^^ (^/- dH)) ^u'^s with a/L in the following way 

9l. = 9:i + {C.9:^)"\l)"^ . (47) 

The coupling to the three dimensional modes produces the usual evolution of in SU{2) 
QCD, where the coupling g'^{a) at the lattice scale a is related to the coupling gQ defined at 
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the infrared scale L as follows 



9\a) = /l . (48) 



Combining the eqs. (^T]) and ( ^8|) we can synchronize the approach to the light cone, i.e. 
the limit 77 ^ with the continuum limit a ^ 0. The condition that the three-dimensional 
evolution of has to be compatible with the two-dimensional evolution of gl^ towards g^^ 
yields that for a ^ the light cone parameter rj approaches zero as 

, , g*J 22 a ^ L 

This relation is similar to the naive scaling result deduced in section 1 from physics 
arguments besides logarithmic modification. Although at the outset we had as parameters 
Ni = L/a and r] in the Hamiltonian we reduce this multiparameter problem to a problem 
with one single coupling constant which can be chosen as g. 

We can use the energy and the critical behavior of the Euclidean 2 + 1 dimensional system 
to obtain its contribution to the ground state energy of the system in 3 + 1 dimensions. In 
our calculations the energy of the zero mode Hamiltonian (h) near the critical coupling 
is proportional to N]_/rj. Using the scaling of 7] determined from the compatibility of 2- 
dimensional and 3-dimensional dynamics, we can eliminate rj in favor of Nl = L/a and get 
a zero mode energy which grows linearly with the 3-dimensional volume. That means the 
zero mode dynamics becomes relevant for the full problem, i.e. 

{ha) ~ NlNI . (50) 

In reality the zero mode system is coupled to the 3 + 1 dimensional system the extension 
of which is large. So we expect that the couplings of the 2 + 1 dimensional system can 
get modified. Because of the universal dynamics near the infrared fixed point this will not 



change the qualitative form of the volume dependence obtained in eq. (50). The behavior 



of the zero mode Hamiltonian has to be taken into account when we try to solve for the 



infinite momentum frame solutions of the complete Hamiltonian. The discussion [g,|16| until 
now often centers on the order of limits. In order to guarantee a meaningful limit of the 
near light cone theory, one first has to take L —>■ 00 and then one can take 77 — *• 0. This 
is a workable procedure in solvable 1 + 1 dimensional models. In QCD in 3 + 1 dimensions 
it hardly can serve as a prescription. Numerical solutions have to be obtained with a finite 
cutoff and the continuum limit can only be reached approximately via scaling relations. It is 
in this case that the role of zero mass excitations will crucially enter the physics. As one sees 
above the synchronization of the continuum limit with the light cone limit leads to results 
which are independent of the near light cone parameter 77. The independence of the ground 
state energy near the critical coupling gives us back Lorentz invariance which states that 
the physics should not depend on the reference frame. The vacuum energy is not allowed 
to depend on rj. There are dynamical zero modes in near light cone QCD which contain the 
physics of the nonperturbative QCD vacuum, the physics of the gluon and quark condensates. 
The nontrivial structure of the nonabelian gauge theory enters in our calculation through the 
kinetic operator of the rotators which at finite resolution have a finite mass gap vanishing 
only in the continuum limit. 
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The near light cone description of QCD based on a modified axial gauge d-A^ = yields 
a very natural formulation of high energy scattering. Polyakov variables near the light cone 
have already entered the calculations of high energy cross sections in two different ways: For 
the geometrical size of cross sections in ref. the correlations between Polyakov lines 

along the projectile and target directions are important. For this situation an interpolating 
gauge may be useful. For the energy dependence of cross sections an evolution equation of 



Polyakov line correlators approaching the light cone has been discussed in ref. which 
generalizes the BFKL- equation. This situation is close to the treatment in this paper. Now 
the decisive step is to back up our theoretical work on effective transverse QCD dynamics 
near the light cone with phenomenological consequences for high energy scattering. 



Appendix 

A choice of At is proposed by comparing Monte-Carlo weights Wt and Wg which link nearest 
neighbors in time and space respectively 

Wtiifii,ip2) = (</^i(x,i/,r + l)|e-^^^>2(a;,i/,r))y^J(^i)J(¥;2) , (51) 
W^iif^^if,) = {Mx + l,y,T)\e-^^^'''\M^,y,T)) . (52) 




Figure 3: Weights Wt{ipi,ip2) of time-like and Wt{ipi,(p2) of space-like links for different 
At. Figures (a,b,c) correspond to At = a/10, a/2, a. 

The probability distribution for the update of the field (p{x, y, r) is given by the product of 
weights for all links connecting it to neighbors in time and space. So the updating procedure 
for if is more effective when Wt and Wg have similar distributions near the angle cp = 7i/2 
where the Jacobian has its maximum. For At = a/10 {cf. Fig. |^) the spatial weights 
Ws have a much wider distribution than the timelike weights Wt, whereas for At = a {cf. 
Fig. 1^) the situation is reverse. For the intermediate case At = a/2 {cf. Fig. ^d) the widths 
of the distributions are very similar, therefore the updating procedure is optimal. 
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